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Abstract 

We present a numerical study of critical phenomena (including the 
Lue- Weinberg phenomenon) arising for gravitating monopoles in a global 
monopole spacetime. The equations of this model have been recently stud- 
ied by Spinelly et al. in a domain of parameter space away from the critical 
points. 
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I. INTRODUCTION 



Topological defects |3J arise in spontaneously broken field theories with non-trivial topo- 
logical structure of the vacuum. If the symmetry broken is global, the defects have infinite 
energy and thus do not posses a particle-like behaviour. One way to get rid of this problem 
is to introduce gauge symmetries, the most popular example being the 't Hooft-Polyakov 
magnetic monopole @. Although initially constructed within an SU(2) gauge field theory, 
monopoles are predicted by all Grand Unified Theories (GUTs) which contain an unbroken 
U(l) symmetry. Since we observe this unbroken £7(1) symmetry in the universe, this leads 
to the so-called "monopole-problem" , the most popular solution to which is the scenario of 
inflation ||. 

The coupling of field theories to gravity enriches their classical spectrum. For example, 
the Lagrangian describing a triplet of self-interacting scalar fields invariant under a global 
0(3) transformation has no finite energy solutions. However, the incorporation of gravity 
leads to global monopoles j|,|5| which have particle-like structure. 

The most striking feature of classical field theory solutions coupled to gravity is the 
pattern of bifurcations found e.g. for the gravitating monopoles merging into black hole 
solutions 0.0]. 

Recently, a model involving both a SO (3) triplet of Higgs fields as well as an 0(3) triplet 
of Goldstone field was considered in [Q. Coupling the Lagrangian of this model to gravity, 
the authors were able to construct finite mass solutions incorporating both the gravitating 
monopole of and the global monopole of 0]. Because many features of these solutions were 
left open, we reconsider here the classical equations and put the emphasis on several unsolved 
questions, namely : (i) how does the topological defect emerge from the purely gravitating 
magnetic monopole, (ii) what is the domain of existence of the solutions in the space of the 
parameters, (iii) do the solutions bifurcate into black holes solutions, (iv) do these features 
persist for large values of the self-interacting coupling constant. These questions are worth 
studying because it is known that the global monopole has a much stronger gravitational 
field at large distances as compared to that of the local monopole. The reason for this is 
that the space-time of the global monopole is not asymtotically flat, while that of the gauged 
monopole is. However, at short distances the gravitational effects of the global monopole 
are weak as long as its mass is much smaller than the Planck mass [0. On the other hand, 
it is well known that local monopoles stop to exist when their radius becomes comparable 
to the corresponding Schwarzschild radius, which implies that the mass of the monopoles 
is of order of the Planck mass Since we are studying critical phenomena of these 

monpoles, the mass of both the local and the global monopole are of the order of magnitude 
of the Planck mass and the argument that the short distance gravitational field of the global 
monopole is weak doesn't hold anylonger. One might thus expect significant changes of the 
critical behaviour of the local monopole in the spacetime of a global one. 

In order to answer these questions, we consider the model of [Q for generic values of 
the expectation values of the local Higgs fields and of the global Goldstone fields. The 
Lagrangian, the spherically symmetric Ansatz and the classical equations are presented in 
Sect. II. The numerical results are discussed in Sect. Ill and illustrated by means of five 
figures. We give our conclusions in Sec. IV. 
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II. THE MODEL 



A. The Lagrangian 

We consider the following Lagrangian density || : 

£m = ~\f; v F^ - l -D^ a D^ a - \d^ a d^ a - V(<j> a i O (i) 
with the potential 

v(r,e) = \m a - v 2 f + ^rr - v 2 ) 2 , (2) 

the covariant derivative of the Higgs triplet a , a = 1, 2, 3 

D,<p a = d,r - ee ahc A h ^ c (3) 

and the field strength tensor 

F; v = d,A a u - d„A° - ee ahc A\A c v . (4) 

This model contains two triplets of scalar fields: the Higgs fields cj) a with vacuum expectation 
value v and self-coupling A and the Goldstone fields £ a , (a = 1, 2, 3) with vacuum expectation 
value i] and self-coupling X g . It is invariant under a SO(3)i oca i x 0(3) g i o b a i transformation : 

(F)'(x) = R ab {x)ct> a {x) , (C)'(x) = R ab C(x) (5) 

where the "i?-part" of the symmetry is gauged by means of the 5*0(3) Yang-Mills field A a ^ 
with gauge coupling e, while the "i?-part" denotes a global transformation. Note that for 
the Lagrangian (|l|) the two scalar triplets are decoupled. The interaction between these 
fields will be carried out through the coupling to gravity. We thus consider the following 
action : 

where G is Newton's constant, g denotes the determinant of the metric tensor g^ v and 
R — g^Rfiv is the Ricci scalar. In this paper, we are carrying out a detailed numerical 
study of the classical, spherically symmetric solutions of the system of equations which arise 
from the variation of @. 

B. Spherically symmetric Ansatz 

For the metric, the spherically symmetric Ansatz in Schwarzschild-like coordinates reads 

ds 2 = g^dx^dx" = -A 2 {r)N{r)dt 2 + N-\r)dr 2 + r 2 {d8 2 + sin 2 6d^ 2 ) , (7) 
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For the gauge, Higgs and Goldstone fields, we use the purely magnetic hedgehog Ansatz 
HI : 



A r a = A t a = 0, (8) 



1 — u(r) „ . „ 1 — u(r) . „ „ 
A e a = — e/ , A/ = ^ sin 6e 6 a , (9) 



(j) a = v h(r)e r a , C = vf(r)e r a . (10) 
We introduce the following dimensionless variable and coupling constants : 

x = evr , a 2 = 4nGv 2 , (3 2 = A , 0* = hi , q =H (H) 

e z y e z v 

The ratio between the radius of the local monopole core r% oc (ev ) _1 and the global monopole 
core r g oc (yAg^) -1 can be given in terms of these quantities : 

— oc = f3 g q 12 

r g ev 

Note that the notation used here corresponds to the one used in || and differs from the one 
in ||. To avoid confusion, we stress here the crucial differences : 

@: a 2 = l-87rG V 2 , g rr = A , g tt = B 
^:a 2 =4nGv 2 , g rr = N- 1 , g tt = A 2 N 



C. Classical field equations 

Varying (|6|) with respect to the metric fields gives the Einstein equations which can be 
combined to give two first order differential equations for A and /x: 

A' = a 2 Ax(- 2 (u') 2 + (h') 2 + (f') 2 ) (13) 

fi' = a 2 x 2 {{U - < ^)+NK) (14) 
x 2 

with the abbreviations 

1C= \(-M) 2 + W) 2 + (J') 2 ) , (15) 

u= ^^-^ + '4 + S + f (a2 - 1)2 + f (/2 - 92)2 • < i6) 

and N and \i are related as follows: 
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N(x) = 1 - 2a 2 q 2 - 2^ . (17) 

Variation with respect to the matter fields yields the Euler-Lagrange equations, which for 
our model are a set of three second order differential equations: 

(ANu'Y = A( U ^ 2 ~^ + h 2 u) , (18) 
x 2 

(19) 

(x 2 ANh')'= A(2hu 2 + (3 2 x 2 h(h 2 - 1)) , (20) 

(21) 

(x 2 ANf )'= A(2f + p 2 x 2 f(f 2 - q 2 )) ■ (22) 

The prime denotes the derivative with respect to x. In order to solve the system of equations 
uniquely, we have to introduce 8 boundary conditions, which we choose to be : 

^(0) = , u(0) = l , /i(0) = , /(0)=0 (23) 



A(oo) = , w(oo) = , h{po) = \ , f(oo)=q (24) 

Note that close to the origin, ji{x — > 0) ~ — a 2 q 2 x such that N(x — » 0) — > 1. Note further 
that while the decay of the Higgs field function h for x — > oo depends on /3, the decay of the 
Goldstone field function / doesn't depend on (3 g : 

(h - 1) w exp(-v / 2 / 9x) , (f - q) ~ c/x 2 for x -> oo . (25) 

The dimensionless mass of the solution is determined by the asymptotic value /i(oo) = /x^ 
of the function fj,(x) and is given by /ioo/a 2 . 



III. NUMERICAL RESULTS 

A. Gravitating monopoles: q = 

For g = 0, no global symmetry breaking takes place and equation (|22|) is trivially solved 
by / = 0. In this limit, the remaining equations are those of the gravitating magnetic 
monopole studied in 0. For completness, we briefly recall the main properties of these 
solutions. The a = limit corresponds to flat space with N(x) = A(x) = 1 and thus the 
't Hooft- Polyakov monopole is recovered. For non- vanishing a 2 , space-time is curved 
and the gravitating monopole arises smoothly from its flat space limit. While the mass 
of the gravitating monopole given (in our rescaled coordinates ) by /ioo/a 2 decreases as 
expected, the ratio /loo/a increases gradually from zero to one. In particular, the function 
N(x) develops a minimum N m at x — x m which becomes deeper for increasing a. At a 
/3-dependent critical value of a (e.g. a cr (j3 = 0) ~ 1.385, a cr {(3 — 1) ~ 1.145) a degenerate 
horizon forms with N m —> and the non-abelian solution ceases to exist. It bifurcates into 
an embedded extremal Reissner-Nordstrom (RN) black hole solution with magnetic charge 
P = 1 and horizon Xh '■ 
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2/ioo a 2 

N RN {x) = 1 — + — , A RN (x) = 1 , fi oa = x h = a (26) 

X x z 

For < (3 < 0.757, the solutions exist up to a maximal value a max with N m (a max ) ^ 
and reach their critical solution on a second branch of solutions with a cr < a max . This 
second branch disappears for (3 > 0.757 and a cr = a max . For intermediate values of /3, i.e. 
for j3 > Ptr — 7.15, however, the pattern of reaching the critical solution changes |10|,ITT|. 



Lue and Weinberg JT0[ observed that for large enough (3 and increasing a, a second local 
minimum of N(x) starts to form which is located between the origin and the minimum which 
corresponds to the Reissner-Nordstrom horizon. Eventually, the inner minimum dips down 
much quicker than the outer one, such that the critical solution is an extremal, non-abelian 
black hole with mass less than that of the corresponding extremal RN solution. For (3 close 
to fa,., the outer minimum is already quite pronounced at the moment the inner minimum 
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B. Gravitating monopoles: q ^ 

One of the main goals of this paper is to study how the critical phenomena described in 
the previous section change in the presence of a global monopole, i.e. for q > 0. 

This is illustrated in FIGs. 1 and 2 for (3 — (3 g = 1 and a = 0.6. FIG. 1 demonstrates the 
evolution of the mass function fi(x) for different values of the parameter q. fi(x) develops 
a local, negative valued minimum and its asymptotic value decreases with q, in particular 
/Zoo becomes negative for q > 0.7. The evolution of the function N(x) is presented in FIG. 
2. The difference N(oo) — N m , where N m denotes the local minimum of N(x), decreases for 
increasing q and the function becomes monotonically decreasing for q > 0.8. The case of 
equal vacuum expectation values q = 1.0 was studied in ||. Note that the solutions cease 
to exist for q > l/{y/2a) w 1.178. 

The natural question to raise now is in which domain of the a, (3, (3 g , q hyperspace 
monopole solutions exist. For the moment, we limit our analysis to the domain of existence 
in the a — q— plane for fixed (3, (3 g . 

First, we analyse the evolution of the solutions for fixed small values of q < qtr{f3,(3 g ) 
and increasing a. We find a very similar picture as in the q = case : the function N(x) 
develops a minimum which approaches zero for a ~ a cr (q, (3, (3 g ). Our numerical analysis 
suggests that the critical value of a increases for increasing q |Jl2fl , e.g. for (3 — (3 g — 1.0 we 
find : 

acr(0.0, 1.0, 1.0) « 1.145 , a cr (0.2, 1.0, 1.0) w 1.145 , 0^(0.6, 1.0, 1.0) « 1.175 (27) 

We further find that the critical value of a depends only little on /3 g , e.g. for q = 0.1 and 
(3 = 1.0 we obtain : 

a cr (0.1, 1.0, 1.0) « 1.145 , a,,. (0.1, 1.0, 5.0) w 1.146 , a cr (0.1, 1.0, 10.0) « 1.147 (28) 

For small values of (3, we find in analogy to that the solutions exist up to a maximal value 
of a, a max (q, (3, (3 g ) with N m ^ and that from there a second branch of solutions exists up 
to a cr < a max with N m reaching zero at a = a cr . Both, a max and a cr increase with q, e.g. 
for (3 = 0.0, (3 g = 1.0, the values are : 
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a cr (0.0, 0.0, 1.0) « 1.385 , a cr (0.3, 0.0, 1.0) « 1.399 , a cr (0.4, 0.0, 1.0) « 1.434 (29) 

For g > q tr (/3,f3 g ) the scenario is quite different as can be guessed from FIG. 2. Indeed, 
when the expectation value of the global Goldstone field becomes large, the function N(x) 
decreases monotonically from N(0) = 1 to its asymptotic value iV(oo) = 1 — 2a 2 q 2 . No 
local minimum develops and the solution just stops existing because the asymptotic value 
iV(oo) = 1 — 2q 2 a 2 of N(x) itself becomes negative. The domain of existence of solutions in 
the a — q— plane is presented in FIG. 3. This FIG. suggests clearly that qt r (P, Pg) decreases 
with decreasing /3, which is indeed what our numerical simulations confirm. 

In order to understand the pattern of reaching the critical solution, we show in FIG. 4. 
the critical solution for q = 0.2. Our numerical results suggest that a degenerate horizon 
forms at Xh ~ 1.148. For > x > Xh, the Goldstone field function f(x) vanishes, while 
all other functions are non-trivial. For x > Xh in contrast, u(x) = 0, h(x) = 1, while the 
Goldstone field function f(x) and the metric functions N(x) and A(x) remain non-trivial 
in this region. We can make a rough approximation to explain this result analytically 
as follows : at and just outside the horizon, we can treat the Goldstone field as roughly 
vanishing f(x) ~ 0. Thus equation ( |T4D becomes 

1 q 2 

H' = a V(— - — + (3 g q A ) for x » x h (30) 
For q = the solution is clearly the RN solution. For q ^ 0, however, we find using (|17f) : 

N(x) = 1 + ^ + - - \a 2 q 4 p 2 g x 2 (31) 

X X o 

where C is an integration constant. This solution has degenerate horizons with N(xh) = 
N'(x h ) = at 



^ 2A4 ' = ± 2^%f 7 ^^ < 32 » 

For our choice of parameters in FIG. 4, one of these four horizons is Xh ~ 1.148. This is 
exactly equal to the horizon we find in our numerical calculations. This approximation is, of 
course, only valid close to the horizon. Away from the horizon, the solution is non-abelian. 
We thus observe a "black hole inside a global monopole". Outside the core of the global 
monopole, where / has reached its vacuum expectation value f(x) = q, (|T4]) reduces to 



the equation for the RN case and using (|17D we obtain for N(x) the metric function of an 
extremal RN black hole with deficit angle 1 — 2a 2 q 2 : 

N RN {x) = l-2a 2 q 2 - 2 -^ + - 2 . (33) 

X X 

The mass of this solution is given by Hoo/a 2 = a~ l \/l — 2a 2 q' 2 . For our choice of parameters 
in FIG. 4, our numerical results indicate that indeed at the critical value of a = 1.145, this 
mass is obtained. 

Finally, we demonstrate in FIG. 5 that the so-called Lue- Weinberg (LW) phenomenon 
observed previously only in the gravitating monopole case for large enough values of j3 
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flO|jrT|| persists in the presence of a global monopole. We show the metric function N(x) for 
(3 = 15, (3 g = 1.0 and q = 0.5. Apparently, with increasing a, first a RN type horizon forms 
at roughly (using (|32|)) x ~ 0.82. At a = 0.7975, a second minimum starts to develop, which 
dips down much quicker than the RN type minimum. At a cr ~ 0.79895, this inner minimum 
has reached zero, while the outer one is still greater than zero. Thus, an extremal non-abelian 
black hole has formed. This scenario is similar to the one observed in the model without 
a global Goldstone field. Comparing the critical value of a, we find that a cr (q,/3 = 15.) 
decreases slightly with increasing q. We find a cr (0.0) = 0.80017 and a cr (0.5) = 0.79895. 



IV. CONCLUSIONS 

We have studied gravitating magnetic monopoles in the spacetime of a global monopole 
and have put emphasis on the study of critical phenomena. We find that the solutions merge 
into extremal black hole solutions representing "black holes inside a global monopole" for 
a choice of parameters where the radius of the core of the local monopole is smaller than 
that of the global one (i.e. for small values of the product f3 g q). The critical value of a 
is increasing with increasing q. For intermediate values of the Higgs boson mass, we find 
that the Lue- Weinberg phenomenon persists in the presence of a global monopole with the 
critical value of a decreasing with increasing q for this phenomenon. 

It is remarkable that many properties of the gravitating monopole persist in the presence 
of a global monopole. Of course, the model studied here involves many parameters and 
we limited our analysis to some particular cases. We believe though that all qualitative 
properties of the solutions are exhibited in our results. 

It would be interesting to study the corresponding black hole solutions of this model, es- 
pecially to investigate how the domain of existence in the a-x^-plane chances in the presence 
of a global monopole. 



Axially symmetric SU(2) monopoles in curved space have been studied in It was 
found that in contrast to flat space, an attractive phase can exist for specific choices of the 
coupling constants. It is left as future work to study the influence of the global monopole 
on the attraction between like-charged gauged monopoles. 
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FIG. 1. The mass function /j,(x) is presented for (3 = (3 g = 1.0, a = 0.6 and for several values 
of q as a function of the dimensionless coordinate x = evr. 
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FIG. 2. Same as FIG. 1 for the metric function N(x). 
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FIG. 3. The domain of existence of the solutions in the a — q— plane is shown for (3 g = (3 = 1.0 
(dashed) and j3 g = 1.0, (3 = (a max : dotted-dashed, a cr : dotted). The solid line represents 
1 — 2q 2 a 2 = 0. The solutions exist below the solid line and to the left of the corresponding dashed, 
dotted-dashed, dotted line, respectively. 
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FIG. 4. The functions N, A, u, h, f are presented for (3 = g = 1.0, q = 0.2 and a w a cr w 1.145 
as functions of the dimensionless coordinate x = evr. 
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